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ABSTRACT 

We have obtained high and intermediate resolution optical spectra of the black-hole 
candidate Nova Sco 1994 in May/ June 1998, when the source was in complete (X-ray) 
quiescence. We measure the radial velocity curve of the secondary star and obtain 
a semi-amplitude of i(r2=215.5±2.4kms^^, which is 6 per cent lower than the only 
previously determined value. This new value for K2 thus reduces the binary mass 
function to f{M) = 2.73±0.09 Mq. Using only the high resolution spectra we constrain 
the rotational broadening of the secondary star, usin i, to lie in the range 82.9- 
94.9 kms~^(95 per cent confidence) and thus constrain the binary mass ratio to lie 
in the range 0.337-0.436 (95 per cent confidence). We can also combine our results 
with published limits for the binary inclination to constrain the mass of the compact 
object and secondary star to the ranges 5.5 - 7.9 and 1.7 - 3.3 MQrespectively (95 per 
cent confidence). Finally, we report on the detection of the Lithium resonance line at 
6707. 8A, with an equivalent width of 55±8 mA. 
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The soft X-ray transient. Nova Sco 1994 (=GRO J1655-40) 
has been studied extensively over the past 3 years in X- 
rays and at optical wavelengths [Tingay et al. 1995, Bailyn 
et al. 1995a, Harmon et al., 1995, Hjellming & Rupen 1995, 
Bailyn et al. 1995b, Paciesas et al. 1996, van der Hooft et al., 
1997, van der Hooft et al., 1998 (hereafter vdH98)]. Strong 
evidence that the compact object in GRO J1655-40 is a black 
hole was presented by Bailyn et al. (1995b) who initially 
established a spectroscopic period of 2.601 ±0.027 days, and 
determined a mass function f(M) = 3.16 ± 0.15 Mq. An 
improved value of f{M) = 3.24 ± 0.09 Mq was presented 
by Orosz & Bailyn (1997; hereafter OB), who measured the 
radial velocity semi-amplitude of the secondary star to be 
if 2 = 228.2 ±2.2 km s"^ and classified it as an F2-F6iv type 
star. Both vdH98 and OB analyzed their quiescent optical 
light curves of Nova Sco 1994 which they combined with the 
OB value for K2 to show that the black hole has a mass 
(Ml) in the range 6.29-7.60 Mq. 

However, in calculating the radial velocity semi- 
amplitude, OB used both quiescent data (taken in 1996 
February 24-25) and outburst data (taken in 1995 April 
30-May 4). Using outburst data in this way can lead to an 
incorrect result as substantial X-ray heating of the secondary 



Table 1. NTT log of Nova Sco 1994 observations 

Date Exp. time Resolution Wavelength range 

29/5/1998 12 X 1800s 0.83 A 6131-6764 A 

30/5/1998 2 X 1800s 4.2 A 4452-7015 A 

31/5/1998 2 X 1800s 4.2 A 4452-7015 A 

2/6/1998 13xl800s 0.83 A 6131-6764 A 



star shifts the 'effective centre' of mass of the star (see Wade 
& Home 1988 and Phillips, Shahbaz & Podsiadlowski 1998). 
This results in a significant distortion of the radial velocity 
curve and hence a spuriously high velocity semi-amplitude. 
The masses of the binary components derived from this will 
therefore be incorrect. 

In this paper we determine the radial velocity curve of 
the secondary star in Nova Sco 1994 using only data taken 
when the system was in quiescence. From our high resolu- 
tion spectroscopy we determine the rotational broadening 
of the secondary star, which, when combined with the new 
value for the mass function and limits to the binary incli- 
nation (vdH98). allows us to determine more precisely the 
individual component masses. 
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Table 2. Sine-wave fits to the radial velocity curvs 



Table 3. Results for the optimal subtraction of the companion 
star 
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Taken from van der Hooft et al. (1998) and fixed. 











1-cr uncertainties are given, see text. 

2 OBSERVATIONS AND DATA REDUCTION 

Intermediate and high resolution optical spectra of Nova Sco 

1994 were obtained on 1998 May 28, 29, 30 and June 1 with 
the 3.5-m New Technology Telescope (NTT) at the Euro- 
pean Southern Observatory (ESO) in Chile using the ESO 
Multi Mode Instrument (EMMI). We used the red arm with 
an order-separating OG 530 filter and gratings #6 and #8 
which gave dispersions of 0.31 A per pixel and 1.26 A per 
pixel respectively (sec Table 1). The TEK 2048x2048 CCD 
was used, binned with a factor two in the spatial direction 
in order to reduce the readout noise. The dispersion direc- 
tion was not binned. Very good seeing allowed us to use a 
slit width of O'.'S which resulted in spectral resolutions of 
0.83 Aand 4.2 A for the gratings #6 and #8 respectively 
(see Table 1). Cu-Ar arc spectra were taken for wavelength 
calibration. Template field stars of a variety of spectral types 
were also observed with intrinsic rotational broadening much 
less than the resolution of our high resolution data. 

The data reduction and analysis was performed us- 
ing the Starlink figaro package, the Pamela routines of 
K. Horne and the molly package of T. R. Marsh. Removal 
of the individual bias signal was achieved through subtrac- 
tion of the mean overscan level on each frame. Small scale 
pixel-to-pixel sensitivity variations were removed with a flat- 
field frame prepared from observations of a tungsten lamp. 
One-dimensional spectra were extracted using the optimal- 
extraction algorithm of Horne (1986), and calibration of the 
wavelength scale was achieved using 5th order polynomial 
fits which gave an rms scatter of 0.02 A. The stability of 
the final calibration was verified with the OH sky line at 
6562. sA whose position was accurate to within 0.1 A. 



3 THE RADIAL VELOCITY OF THE 
COMPANION STAR 

Our first task was to measure the radial velocity of the 
F-type secondary star. The radial velocities were mea- 
sured from the spectra by the method of cross-correlation 
(Tonry & Davis 1979) with a template star. Prior to cross- 
correlation the spectra were interpolated onto a logarithmic 
wavelength scale (pixel size 14.5 kms~^) using a sin x/x in- 
terpolation scheme to minimize data smoothing (Stover et 
al. 1980), and then normalised. We masked the interstellar 
feature at 6280 A and Ha before correlating the spectra. We 
used HR2927 as our template star after we had broadened it 
by 90 kms~^in order to simulate the spectrum of Nova Sco 



1994 (a more accurate value for the rotational broadening 
is derived in the next section). The radial velocity of the 
template star (derived using the position of the Hq absorp- 
tion line to be —6.7 kms^^) was then added to the radial 
velocities of Nova Sco 1994. 

Using the orbital ephemeris given by vdH98 we phase- 
folded the heliocentric radial velocities and then fitted a sine 
wave (see Table 2 and Figure 1), obtaining a of 0.99. We 
then combined the data with the quiescent radial velocity 
data of OB and repeated the fitting procedure obtaining 
a significantly worse value of xS=4.9. Note that our high 
resolution data has a factor of 4 better spectral resolution 
than that of OB; the uncertainty in our individual radial 
velocity measurements being typically 6kms~^. 



4 THE SPECTRAL TYPE AND ROTATIONAL 
BROADENING OF THE COMPANION STAR 

We determine the spectral type of the companion star by 
minimizing the residuals after subtracting diff'erent template 
star spectra from the Doppler-corrected average spectrum. 
This method is sensitive to the rotational broadening v sin i 
and the fractional contribution of the companion star to the 
total flux. The template stars we use are in the spectral 
range F0-F8, V-III and were obtained during this observ- 
ing run but also from previous runs at La Palma and with 
comparable dispersion. 

First we determined the velocity shift of the individual 
high resolution spectra of Nova Sco 1994 with respect to 
each template star spectrum, then Doppler-averaged to the 
rest frame of the template star (see section 3). In order to 
determine the rotational broadening, wsin i, we follow the 
procedure described by Marsh, Robinson & Wood (1994). 
Basically we subtracted a constant, /, representing the frac- 
tion of light from the template star, multiplied by a rotation- 
ally broadened version of that template star. The optimal 
values of v sin i and / are obtained by minimising (see 
Table 3). The above analysis was performed in the spectral 
ranges 6380-6520 A and 6600-6750 A which excludes Ha. 
This was the only region common to all the templates stars 
and the high resolution spectra of Nova Sco 1994. We used 
a linear limb-darkening coefficient of 0.52 (Al-Naimiy 1978) 
appropriate for 6500 A and an effective temperature of 6500 
K (typical for a F star) . Figure 3 shows the results for those 
templates with / < 1. It should be noted that the main se- 
quence template stars all gave / > 1, as also found by OB. 
From Table 3 it can be seen that the minimum occurs 
at spectral type F6 with a v sin i of 91.2±3.7 kms"^ The 
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Table 4. The rotational broadening of the companion star and 
the binary mass ratio. 90 per cent confidence levels are given. 

Limb darkening vsin i q{=M2/Mi) 

coefficient (u) (km s~^) 

0.00 86.5 (+3.3, -3.6) 0.360 (+0.028, -0.033) 

0.52 91.2 (+3.7, -3.7) 0.400 (+0.036, -0.030) 

90 per cent confidence levels shown were obtained using a 
Monte Carlo simulation (see section 5 for details), after the 
x1 of the fit was multiplied by a scaling factor so that the 
X?,min was 1 (Lampton et al. 1976). 

The analysis above assumes that the limb-darkening co- 
efficient appropriate for the radiation in the line is the same 
as for the continuum. However, in reality this is not the 
case; the absorption lines in early-typo stars will have core 
limb-darkening coefficients much less than that appropri- 
ate for the continuum (Collins & Truax 1995). In order to 
determine the extreme limits for vsin i we also repeated 
the above analysis for the F6 template star using zero limb- 
darkening and found u sin 1=86^3 ^ km s~^ (90 per cent con- 
fidence). Wo thus constrain wsin i to lie in the range 82.9- 
94.9 kms"^(95 per cent confidence; since we have two one- 
tailed distributions representing the upper and lower limits). 

As can be seen from Figure 3 (top) the determination 
of the spectral typo depends critically on the choice of tem- 
plate spectral type and luminosity class. Even if we had 
many template stars we would find that using stars of the 
nominally same spectral class and type would give differing 
results. In Figure 3 we also show a fit to the data; this is 
what would be expected if there were no systematic effects. 

In order to estimate this systematic effect in our analysis 
of Nova Sco 1994 we performed a Monte Carlo simulation 
involving the same analysis but now using a a template star 
of known spectral type (F5iv; HR870) as the target. To this 
target we added noise to produce a spectrum of comparable 
quality to our Doppler-averaged Nova Sco 1994 spectrum. 
We then repeated the broadening and optimal subtraction 
procedure using the same 21 templates star as was used for 
the Nova Sco 1994 analysis, thereby determining the best fit. 
We found that we could determine the spectral typo of the 
trial template star to within 1 sub-class. Figure 3 (bottom) 
shows a single simulation; note the systematic effects, i.e. 
the scatter within each spectral type bin. 



5 THE BINARY SYSTEM PARAMETERS 

Since the companion star fills its Roche lobe and is syn- 
chronised with the binary motion, the rotational broaden- 
ing provides a direct measurement of the binary mass ra- 
tio, q (=M2/Mi), through the expression usin i = [-fi'2(l + 
g)0.49g2/=V[0-6g^/^ +ln(l + 5)1/3] (Home, Wade & Szkody 
1986). Substituting our values for K2 and v sin i we calculate 
q (see Table 4). In order to determine the uncertainty in q we 
used a Monte Carlo simulation, in which we draw random 
values for the observed quantities which follow a given dis- 
tribution, with mean and variance the same as the observed 
values. For K2 the random distribution to taken to be Gaus- 
sian as the uncertainty is symmetric about the mean value. 



Table 5. The masses of the binary components. 90 per cent con- 
fidence levels are given. 
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0.52 
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However for v sin t the uncertainty is asymmetric and so we 
had to determine the actual distribution numerically. This 
was done by first calculating the maximum likelihood dis- 
tribution using the actual fit values (after rescaling the 
xl values so that the xl,min was 1) and then determining 
the cumulative probability distribution. By picking random 
values (from a uniform distribution) for the probability we 
obtain random values for vsin i. Given the uncertainty in 
the limb darkening coefficient and hence vsin i, we find q 
to lie in the range 0.337 - 0.436 (95 per cent confidence). 

Using our values for K2 and q with the orbital period, P, 
and the binary inclination, i, wc can determine the masses 
of the compact object, Afi and the companion, M2 star us- 
ing PKI/2ttG = Ml sin^ i/(l + g)^. The binary inclination 
is given by vdH98, and lies in the range 63.7° and 70.7°. 
It should be noted that the solutions for the binary in- 
clination and hence binary masses, determined by vdH98 
depend on the flaring angle of the accretion disc. Here we 
use the range for the binary inclination which encompasses 
the range of flaring disc angles assumed by vdH98. Table 
5 shows the values wc obtain for the masses which depend 
on the limb-darkening coefficient used to determine v sin i 
and the binary inclination. The 90 per cent confidence level 
was determined using a Monte Carlo simulation (sec Figure 
4) similar to that described above. Given the uncertainties 
in the limb-darkening coefficient (and thus in v sin i and q) 
and the binary inclination, we find Mi and M2 to lie in the 
range 5.5 - 7.9 Meand 1.7 - 3.3 M© respectively (95 per cent 
confidence) . 

By modelling the optical light curves of Nova Sco 1994 
vdH98 and OB have also determined q. They find q to lie 
in the range 0.24-0.42 (3-cr) and 0.289-0.385 (3-fr) respec- 
tively. (However, it should be noted that the errors given 
by OB arc only internal statistical errors, whereas the anal- 
ysis of vdH98 includes the systematic uncertainties.) Us- 
ing optical spectroscopy (see section 5), we find q to lie 
in the range 0.337-0.436 (95 per cent confidence). Also 
from optical spectroscopy, the projected radial velocity semi- 
amplitude for the compact object has been determined by 
Soria et al., (1998) to be 76.2±7.5 kms"^ We can use this 
along with our value for K2 (see section 3) to determine 
g; g=0. 354+0.035 (l-cr). As one can see, the four indepen- 
dent measurements for q above are fully consistent with each 
other. This consistency is also seen in the determination of 
the binary inclination; van der Hooft 1997, vdH98, and OB 
determine values for i which are all also consistent with each 
other. 
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6 THE DETECTION OF LITHIUM? 

High lithium abundances appear to be a common feature 
of late-type secondary stars in neutron stars and black hole 
SXTs, a result with is totally unexpected given the post- 
main sequence dilution and the mass transfer history of these 
stars (Martin et al., 1994). In order to counteract the deple- 
tion of lithium, caused by convectivc mixing, it is necessary 
to invoke a process to produce lithium. Several mechanisms 
have been put forward to explain the high lithium abun- 
dance in the neutron star and black hole SXTs; lithium syn- 
thesis in the supernova explosion that formed the compact 
object (Woosley et al., 1990), or a — a reactions during the 
strong SXT outbursts (Martin et al., 1994). 

We can use the ratio of Li l/Ca l equivalent width 
(^RLi/Ca), which is not sensitive to veiling for deriving up- 
per limits to the surface Li abundance (see Figure 5). For 
the other SXTs A0620-00, V404 Cyg and Cen X-4 RLi/ca 
has been found to be > 1 (Martin et al., 1994); for Nova Sco 
1994 we find RLi/ca=0.48±0.08. The equivalent widths of 
Li I and Ca l arc 55±8 mA and 114±10 mA rrespectivoly. 

The interpretation of the lithium abundance in Nova 
Sco 1994 is somewhat difficult. The secondary stars in SXTs 
are tidally locked and therefore have high rotational veloc- 
ities, hence, convectivc depletion may be inhibited by re- 
ducing the angular momentum transport in the base of the 
convection zone. One might expect there to be a correlation 
between rotation rate and lithium abundance (Pallavicini, 
Cerruti-Sola & Duncan 1987), since convective mixing may 
cause additional depletion if lithium is taken to layers where 
it can be burnt. However, the link between rotation and 
lithium in subgiant stars is not yet well established; de 
Medeiros, do Nascimento Jr & Mayor (1997) have found 
that the lithium abundance is independent of rotational ve- 
locity. Indeed the sample tabulated by Balachandran (1990) 
shows a remarkable range of lithium equivalent width (2- 
100 mA) within F6iv stars but there is no correlation with 
rotational velocity. Thus it is difficult to compare the actual 
lithium equivalent width in Nova Sco 1994 with that in "nor- 
mal" F6iv stars, such as our template star HR2927 in Fig. 
5. All we can say is that we have detected lithium in Nova 
Sco 1994, but whether it is intrinsic to the secondary star 
remains unclear. Spallation appears to bo the most obvious 
mechanism as it is the most probable cause of the high Li 
abundance seen in other SXT, but in Nova Sco 1994 mech- 
anisms related to the spectral type, luminosity of the star 
and magnetic braking have to be considered as well. 



7 DISCUSSION 

From Figure 2 it can be seen that after removal of the sec- 
ondary star from the average spectrum of Nova Sco 1994, 
there is narrow residual Ha emission (EW=10.0±0.5A; 
FWHM=146 kms^^) which could be either due to chromo- 
spheric activity on the secondary star, and/or X-ray heat- 
ing on the inner face or emission from the bright spot. Ha 
emission arising from the quiescent accretion discs in SXTs 
usually have a FWHM of a few thousand kms~^; the resid- 
ual feature we see has a FWHM of only a few hundred 
km s"'^ (e.g. Casares et al., 1997), so it must arise from re- 
gions on the secondary star or the outer edge of the accre- 



tion disc (Marsh et al,, 1994). Our data do not have suf- 
ficient phase resolution to allow us to distinguish between 
the above possibilities; only a detailed phase-resolved radial 
velocity study of the narrow component will solve its origin. 

However, it is interesting to note that the narrow omis- 
sion can be powered solely by X-ray heating. An upper limit 
to the X-ray flux on Nova Sco 1994 is set by RXTE ASM 
observations &t < Fx = 3.1 x 10~^° erg cm^ s~^. The X-ray 
irradiation at the secondary star would then be < 2.3 x 10^" 
erg cm s^ , where wc have used a distance of 3.2 kpc 
(Hjellming & Rupen 1995) and a binary separation of 16.6 
R©. At such levels the secondary would be receiving suffi- 
cient energy to power the narrow Hq emission. 

In calculating the radial velocity semi-amplitude of the 
secondary star OB used both quiescent (1996 February 24- 
25) and outburst data (1995 April 30-May 4), and fitted 
the resulting curve with a sine wave. It has been known, 
especially in studies of dwarf novae (Wade & Hornc 1988 and 
Davey & Smith 1992) that using outburst data in this way 
leads to incorrect results. The effect of substantial heating 
of the secondary causes the centre of light, given by the 
strength of the absorption lines, to shift from the centre of 
mass. This results in a significant distortion of the radial 
velocity curve leading to a spuriously high radial velocity 
semi-amplitude (Davey & Smith 1992). The binary masses 
derived from this will therefore be incorrect. 

Phillips, Shahbaz & Podsiadlowski (1998) have esti- 
mated K2 by fitting only the outburst radial velocity data of 
OB. They use a crude model based on X-ray heating of the 
secondary star and estimate K2 to lie in the range 194-214 
kms~^(90 per cent confidence), assuming ^=0.33. Note that 
this range compares well with our value of /f2=215.5±2.4 
km s""*^ obtained by fitting our high resolution radial velocity 
data of Nova Sco 1994 in complete quiescence. We therefore 
conclude that the masses of the Nova Sco 1994 system are 
somewhat lower than the OB value, but the compact object 
mass is still sufficiently high (> SM©) to imply that it is a 
black hole system. 
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Figure 1. Top panel: The radial velocity curve of the secondary star obtain using only our NTT data. The solid line shows a sinusoidal 
fit to the data. Bottom panel: The radial velocity curve of the secondary star obtain using our NTT data (crosses) and also the data of 
OB (stars). The solid lines in each panel show a sinusoidal fit to the data. The data have been folded on the orbital ephemeris given by 
van der Hooft at al., (1998), and 1.5 cycles orbital are shown. 
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Figure 2. The results of the optimal subtraction. From top to bottom: the residual spectrum of Nova Sco 1994 after subtracting the 
template star times /=0.98, the variance-weighted Dopplcr averaged spectrum of Nova Sco 1994, the template F6iv star (HR2927) 
broadened by 91.2 kms~^, and the template F6iv star again to show the narrow absorption lines. The spectra have been normalized and 
shifted vertically for clarity. 
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Figure 3. Top: The xi of optimal subtraction analysis as a function of the spectral type of the template star used. The solid line is a 
quadratic fit to the data. Bottom; Same as above but for a sample star of known spectral type (F5). Different symbols indicate different 
luminosity classes: triangles for dwarfs as shown. Only template stars that gave / < 1 in the analysis are shown (see text). 
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Figure 4. The mass of the binary components obtained using a Monte Carlo sinmlation (see section 4). 90 per cent confidence regions 
are shown. The sohd elhpses are regions obtained using the continuum value for tlic limb-darlicning coefficient, where as the dashed 
ellipses are regions obtained using no limb darkening. The top and bottom two ellipses are for a binary inclinations of 63.7° and 70.7° 
(van der Hooft et al., 1998). 
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Figure 5. From top to bottom; The Doppler averaged spectrum of Nova Sco 1994, a template star (F6iv) broadening by 91 kms ^to 
match Nova Sco 1994, the same template, to show positions of the absorption lines more clearly. The spectra have been offset for clarity. 



© 0000 RAS, MNRAS 000, 000-000 



